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Interactions between swimming cells and surfaces are essential to 
many microbiological processes, from bacterial biofilm formation to 
human fertilization. However, in spite of their fundamental impor- 
tance, relatively little is known about the physical mechanisms that 
govern the scattering of flagellated or ciliated cells from solid sur- 
faces. A more detailed understanding of these interactions promises 
not only new biological insights into structure and dynamics of flag- 
ella and cilia, but may also lead to new microfluidic techniques for 
controlling cell motility and microbial locomotion, with potential ap- 
plications ranging from diagnostic tools to therapeutic protein syn- 
thesis and photosynthetic biofuel production. Due to fundamental 
differences in physiology and swimming strategies, it is an open ques- 
tion whether microfluidic transport and rectification schemes that 
have recently been demonstrated for pusher-type microswimmers 
such as bacteria and sperm cells, can be transferred to puller-type 
algae and other motile eukaryotes, as it is not known whether long- 
range hydrodynamic or short-range mechanical forces dominate the 
surface interactions of these microorganisms. Here, using high-speed 
microscopic imaging, we present direct experimental evidence that 
the surface scattering of both mammalian sperm cells and unicellular 
green algae is primarily governed by direct ciliary contact interac- 
tions. Building on this insight, we predict and verify experimentally 
the existence of optimal microfluidic ratchets that maximize recti- 
fication of initially uniform Chlamydomonas reinhardtii suspensions. 
Since mechano-elastic properties of cilia are conserved across eu- 
karyotic species, we expect that our results apply to a wide range of 
swimming microorganisms. 



Surface interactions of motile cells play crucial roles in a 
wide range of microbiological phenomena, perhaps most 
prominently in the formation of biofllms pQ and during the fer- 
tilization of mammalian ova [2 . Yet, in spite of their widely 
recognized importance, the basic physical mechanisms that 
govern the response of swimming bacteria, algae or sperma- 
tozoa to solid surfaces have remained unclear. This predica- 
ment is exemplified by the current debate [3] SI [5j [6] about 
the relevance of hydrodynamic long-range forces and steric 
short-range interactions for the accumulation of flagellated 
cells at liquid-solid interfaces. From a general perspective, im- 
proving our understanding of cell-surface scattering processes 
promises not only new insights into structure, dynamics and 
biological functions of flagella and cilia - it will also help to 
advance microfluidic techniques for controlling microbial lo- 
comotion [7, 8 , with potential applications in diagnostics [9], 
therapeutic protein synthesis [TO] and photosynthetic biofuel 
production [111 1121 IT3l 114] . That microfluidic circuits provide 
an excellent testbed for developing and assessing new strate- 
gies for the control of cell motility was recently demonstrated 
by the rectification of random bacterial swimming through mi- 
croscopic wedge-shaped barriers [7, 8 . However, since eukary- 
otic and prokaryotic swimming strategies differ substantially 
from each other [6, 15, 16, 17l ll8lfl9] . it is unclear whether de- 
sign principles that exploit surface collisions to achieve control 
of bacterial locomotion are transferrable to motile eukaryotes. 

Aiming to elucidate the role of eukaryotic cilia in cell- 
surface interactions, we report here a detailed experimen- 
tal investigation of surface scattering for bull spermato- 
zoa and Chlamydomonas reinhardtii algae (simply referred 



to as Chlamydomonas in the following). Bull sperm and 
other mammalian spermatozoa are "pusher" swimmers that 
generate propulsion by undulating a single posterior cilium 
(Fig. [IK). By contrast, a wild-type Chlamydomonas cell is 
a "puller" that achieves locomotion by the breaststroke-like 
beating [T5J [161 [17] of a pair of anterior flagella (Fig. |2|\). 
Chlamydomonas algae have long been appreciated as premier 
model organisms in biology [20 , 21 , 22 , in particular for study- 
ing photosynthesis [TTJ [23] and ciliary [15] [20] [24] functions 
in eukaryotes. More recently, they have also attracted consid- 
erable interest as possible sources of therapeutic proteins [10 
and renewable biofuels [HIIIlElIllEllESlllllEH]. Against 
this backdrop, our second goal is to demonstrate the feasibility 
of microfluidic rectification schemes for these organisms. 

Rectification of bacterial run-and-tumble motion in mi- 
crofluidic ratchets 7, 8 is believed to result from the swim- 
mers' tendency to align their motion along the ratchet barri- 
ers, either by steric [4, 5, 6 or by hydrodynamic 3 surface 
interactions, although the exact mechanism is not well under- 
stood. While steric alignment with surfaces seems intuitively 
plausible for rod-shaped bacteria like E. coli or B. subtilis, ad- 
ditional hydrodynamic alignment is thought to arise from the 
fact that the posterior flagellar bundle of a bacterium creates 
a "pusher" -like dipolar flow during locomotion [19]. This flow 
points outwards along the body axis and inwards along the 
lateral directions [6] . The presence of a wall couples the swim- 
mer's translation and rotation, and causes it to align parallel 
to the surface [3]- By contrast, the anterior flagella of wild- 
type Chlamydomonas "pull" the organism through the fluid, 
thereby generating a far- field flow topology [161 [T7] that looks 
roughly opposite to that of a bacterium. Hence, far- field hy- 
drodynamics suggests that Chlamydomonas should either turn 
away from or collide head-on with a nearby no-slip surface, but 
the complex time-dependent flow structure [161117] close to the 
cell body makes it difficult to predict the scattering dynamics 
in the vicinity of the surface. It is therefore not possible to 
infer from general hydrodynamic arguments whether it is at 
all feasible to design microfluidic structures that are capable 
of rectifying algal swimming. Moreover, purely hydrodynamic 
considerations completely neglect direct contact interactions 
between cilia or flagella and solid surfaces. Unfortunately, 
this potentially important scattering mechanism [4] is not in- 
cluded in currently prevailing theoretical models of microbial 
swimming near solid boundaries [3]. 
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Here, we present direct experimental evidence that the 
scattering of bull spermatozoa and Chlamydomonas algae off 
a solid boundary is, in fact, mainly determined by the con- 
tact interactions between their flagella and the surface, while 
hydrodynamic effects only play a secondary role. Building 
on these insights we derive a simple criterion to predict an 
efficient ratchet design for Chlamydomonas and confirm its 
validity experimentally, thereby demonstrating that robust 
rectification of algal locomotion is possible. More generally, 
our results show that the interactions between swimming mi- 
croorganisms and surfaces are more complex than previously 
recognized, suggesting the need for a thorough revision of cur- 
rently accepted paradigms. Since mechano-elastic properties 
of eukaryotic cilia are conserved across eukaryotic species, we 
expect flagella-surface interactions to play a similarly impor- 
tant role for a wide range of natural micros wimmers, thus 
promising new diagnostic tools and microfluidic sorting de- 
vices for sperm [9] and other motile cells. 

Results 

To identify the dynamical details of eukaryotic cell-surface 
interactions, we studied the surface scattering of bull sper- 
matozoa and four different Chlamydomonas algae strains in 
quasi-2D microfluidic channels (height ~ 25/xm), using high- 
speed microscopic imaging (Materials and Methods). 

Scattering of Individual Sperm Cells from Solid Boundaries. 

We analyzed the scattering of individual sperm cells from 
corner-shaped channel boundaries at high and low temper- 
ature (Fig. 1A,B and Movies SI and S2). Direct observation 
reveals that the short-range interaction of the cilium with the 
boundary determines how a spermatozoon swims along a solid 
surface. In a typical scattering event, a sperm cell closely fol- 
lows the boundary until it reaches a corner and departs at an 
angle 6, defined here relative to the initial swimming direction 
such that = corresponds to the surface tangent (Fig. IB). 
We determined the temperature-dependent distributions of 
from more than 200 scattering events by tracking the position 
of the cell body up to a distance of 70 /im from the corner 
(Fig. 1C). The histograms show that sperm-surface interac- 
tions are typically characterized by negative scattering angles 
< 0, due to the fact that the beat amplitude of the ciliary 
motion is much larger than the size of the cell body. Hence, 
ciliary contact with the surface tends to turn the spermatozoa 
towards the boundary, thereby preventing their escape from 
fiat and weakly curved surfaces. 

To illustrate how differences in the ciliary beating pat- 
terns affect the surface interaction of spermatozoa, we exploit 
the fact [29] that amplitude and frequency of the ciliary mo- 
tion increase monotonously with temperature T in the range 
5°C < T < 38 °C. This variability is caused by a change in 



motor activity [30], approximately described by an Arrhenius 
law oc exp(— A/kT) [29], where A denotes the activation en- 
ergy and k the Boltzmann constant. If ciliary contact governs 
the surface interactions of sperm, then one should expect that 
the absolute mean scattering angle \0\ increases with temper- 
ature. By comparing the scattering distributions at low and 
high temperatures, we find that this is indeed the case (Fig. 
1C). At low temperature T = (10 1)°C the experimental 
data yields a mean scattering angle = (— 5.6±1.0)°, whereas 
= (-12.6 ± 0.7)° at high temperature T — (29 ± 1)°C. In 
particular, these results suggest that, at higher temperatures, 
sperm cells can become more easily trapped at strongly curved 
surfaces. With regard to future biotechnological applications, 
this self-trapping by ciliary beating can provide a useful mech- 
anistic basis for sorting and rectifying spermatozoa [9]. 

Flagella-lnduced Scattering of Individual Algae from Solid 
Boundaries. To further test the idea that ciliary contact dom- 
inates eukaryotic cell-surface interactions, we studied the sur- 
face scattering of four different Chlamydomonas strains (The 
Chlamydomonas Resource Center, http : //www . chlamy . org) : 
the wild-type CC-125, the long-flagella mutant CC-2289 lf3- 
2, the short-flagella mutant CC-2347 shfl, and the moving- 
backwards-only mutant CC-2679 mbol (Fig. |2j. All four 
Chlamydomonas strains share an essentially identical geomet- 
rical structure, but while the wild-type and the mutants lf3-2 
and shfl are "puller" -type swimmers that differ only in flagella 
length (6 — 8/xm for shfl 31 , compared to 11 — 13 /im for wild- 
type and 12 — 22 /im for lf3-2 [32]), the backwards-swimmer 
mbol has a persistent undulatory swimming gait and can be 
considered a "pusher" similar to bacteria and spermatozoa 
(see Movies S3, S4, S5, S6, S7). For each strain, we recorded 
more than 1300 boundary-scattering events in quasi-2D mi- 
crofluidic channels (height ~ 25/xm). 

High-speed imaging of individual cell trajectories for both 
wild-type and mutants reveals that the interaction of Chlamy- 
domonas with the channel wall is also strongly affected by 
the short-range contact forces between the flagella and the 
surface (Fig. [2j\,B and Movies S3, S4, S5, S6). In the case 
of the three fuller" swimmers (wild-type, lf3-2 and shfl), 
the flagella prevent the cell body from touching the surface 
while simultaneously creating an effective torque that turns 
the organism away from it (Fig. |2]\,E). By contrast, for the 
backwards-swimming mutant mbol, posterior thrust by the 
flagella pushes the cell body onto the surface (Fig. [2p,E). 
Subsequently, the mechanical contact of the cilia with the 
boundary leads to a net torque that keeps rotating the alga 
towards the surface (Fig. [2fe). As a result, mbol cells remain 
trapped at the channel wall for several seconds, compared to 
< 0.5 s for wild-type. 

The contact force that is exerted by a flagellum onto the 
surface per stroke can be estimated from the flagellar beat 
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Fig. 1. Surface scattering of bull spermatozoa is governed by ciliary contact interactions, as evident from the scattering sequences of individual cells at two temperature 
values (A) T = 10°C and (B) T = 29°C. The background has been subtracted from the micrographs to enhance the visibility of the cilia. The cyan-colored line indicates 
the corner-shaped boundary of the microfluidic channels (see Movies SI and S2 for raw imaging data). The horizontal dotted line in the last image in (B) defines = 0. Scale 
bars 20 /im. (C) The probability distributions of scattering angles 6 from the corner peak at negative angles, due to the fact that the beat amplitude of the cilia exceeds the 
size of the cell body (sample size: n = 116 for T = 10°C and n = 115 for T = 29°C). At higher temperatures, the cilia exhibit a larger oscillation amplitude and beat 
frequency [29], resulting in a larger swimming speed and shifting the typical scattering angles to larger absolute values. 
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frequency (~ 50 Hz) and the experimentally observed angu- 
lar displacement amplitude per beat (~ 0.2 rad), which gives 
a typical angular speed uj ~ 10 rad/s. Assuming a spher- 
ical cell body (radius a ~ 5 /im) , the torque T can be ob- 
tained from T ~ £cj, where the rotational drag coefficient is 
given by £ = 87rrja 3 . Using r\ — 10 -3 Pa-s for water, we 
find T 30pN-/im and, furthermore, by assuming a flag- 
ella length L ~ 10 /im, the typical force F ~ T/L ~ 3pN. 
These estimates, which are based on the observed rotation of 
the cell body near the surface, are consistent with the val- 
ues obtained from recent measurements of freely swimming 
Chlamydomonas algae [171 133] , suggesting that reorientation 
at the wall is primarily determined by flagellar contact. 

To quantify the different surface scattering laws for each 
of the four Chlamydomonas strains in detail, we measured the 
incidence and scattering angles 0i n and ou t (see Materials and 
Methods) and determined the conditional scattering distribu- 
tions P(#out|#m), defined as the probability of being scattered 
into the interval [# ut,#out + d0 out ] for a given incidence an- 



gle Oin (Fig. [2p). For wild- type Chlamydomonas, P(#out|#in) 
is independent of the incidence angle 0i n and exhibits a narrow 
peak at O ut ~ 16°. The typical scattering angle O ut ~ 20° is 
larger for the long-flagella mutant lf3-2, whereas for the short- 
flagella mutant shfl, the maximum of P(# ou t|#in) is shifted to 
smaller angles ou t ~ 12°. The systematic increase of the typ- 
ical scattering angle with flagella length, in conjunction with 
the above force estimates, implies that, in all three cases, the 
characteristic escape angle ou t is set geometrically through 
the length of the cilia, the diameter of the cell body and the 
typical distance of the latter from the surface at the moment 
of departure (see Movies S3, S4, S5). For shfl, the cell body 
can come closer to the wall than for the wild-type, and the re- 
sulting lubrication forces may also affect the escape dynamics. 
The relatively larger spread in the distribution of scattering 
angles for shfl is compatible with previous observations of 
stronger intrinsic fluctuations in shorter flagella |34j . By con- 
trast, the "pusher" -mutant mbol generally remains close to 
the surface, ou t ~ or ou t ~ ty (Movie S6). The differ- 
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Fig. 2. Surface scattering of Chlamydomonas is governed by ciliary contact interactions. (A) Scattering sequence for wild-type Chlamydomonas CC-125 (Movie S3). 
Upper panel: original micrographs. Lower panel: cilia manually marked red. Results for the long-flagella mutant lf3-2 and the short-flagella mutant shfl look qualitatively 
similar (Movies S4 and S5). Scale bar 20 /im. (B) The mutant pusher mbol remains trapped for several seconds (Movie S6). Scale bar 20 /im. (C) The conditional probability 
distributions -P(#out|#in) indicate that, for all four strains, memory of the incidence angle is lost during the collision process, due to multiple flagellar contact with the 
surface. (D) The cumulative scattering distribution P(0 ou t) shows how cilia length and swimming mechanisms determine the effective surface-scattering law. (E) Schematic 
illustration of the flagella-induced scattering and trapping mechanisms. 
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ences between the scattering laws of the four different strains 
are also clearly evident from the mean scattering distributions 
-P(0out), obtained by averaging P(#out|#in) over all incoming 
angles (Fig. |2^1D. The symmetric bimodal shape of P(0 ou t) 
for mbol signals a complete loss information about the inci- 
dence angle, whereas in the case of the three "puller" strains 
the swimmers still remember their incoming directions, but 
not the exact values of 6i n . Quantitatively similar probability 
distributions characterize scattering laws in thicker chambers 
(80 /im and 300 /mi). 



Optimal Rectification of Algal Locomotion in Microfluidic 
Ratchets. Knowledge of flagella-induced scattering can be 
used to design microscopic obstacles that will passively guide 
the random swimming of microorganisms in a desired direc- 
tion (Fig. 3). The design principles for mbol mutants and 
spermatozoa [9] are similar to those for E. coli |8], as 
these species glide along surfaces after collision. The chal- 
lenge is to find the optimal rectification geometry for "puller" 
organisms, which scatter at a finite angle O ut > off solid 
surfaces. To demonstrate the optimization procedure, we fo- 
cus on the wild-type Chlamydomonas strain, which has the 
peak angle # ou t ~ 16°, and consider wedge-shaped obstacles, 
like those successfully employed in bacterial rectification 18] 
(see Fig. 3B). Optimal rectification requires maximizing the 
ratio Jf/jb between the algae currents in forward and back- 
ward directions. The diffusive backward flux jb is driven by 
gradients in the algal concentrations across the obstacle rows 
(Fig. 3C) and can be minimized by decreasing the gap dis- 
tance do between neighboring obstacles relative to the effec- 
tive width of a single barrier, do <C ds (Fig. 3B). At the 
same time, the forward current jp can be maximized by ad- 
justing the wedge angle a of the barriers to exploit secondary 
scattering (Fig. 3B). Assuming a fixed angle O ut 3 basic geo- 
metric considerations yield the criterion 2a + ou t < tt/2 for 
secondary scattering in the forward direction. Thus, for ideal 
deterministic scatterers with ou t #out one would expect 
maximal rectification for a* « 7r/4 — # ou t/2, yielding a* « 37° 
for wild- type Chlamydomonas. This estimate should, how- 
ever, be viewed as an upper bound for the optimal opening an- 
gle, since the measured scattering distributions P(6 ou t) have 
large-angle tails (Fig. [2p) that shift the optimal a- value to 
smaller angles. 

We tested this prediction in numerical simulations of a 
minimal 2D model and by performing quasi-2D microfluidic 
experiments (Fig. 3). In our simulations (see Materials and 
Methods for details), microswimmers are represented by non- 
interacting point particles whose motion captures the main 
observed features of Chlamydomonas trajectories. Each parti- 
cle moves ballistically at a constant speed V, performs random 
turns after exponentially distributed run-periods 15 with 
persistence time r ~ 1.5 s, and scatters from boundaries at 
a random angle out . Both V and ou t are drawn from dis- 
tributions that mimic the experimentally measured speed dis- 
tributions and scattering laws (see Materials and Methods). 
The value of r is consistent with estimates from the mean 
square displacement of wild-type Chlamydomonas, as mea- 
sured in the quasi-2D microfluidic chambers (height 25/im) 
used for the rectification experiments (Fig. 3A). Compared 
with 3D chambers, the typical swimming speed of the algae 
(V ~ 30/xm/s in quasi-2D chambers) is reduced by approxi- 
mately a factor of two due to the presence of nearby no-slip 
surfaces, which also tend to suppress hydrodynamic interac- 
tions [35] . 

In both experiments and simulations, we considered cham- 
bers with four compartments separated by rows of wedge- 



shaped barriers with gaps of length cfc (Fig. 3B,C). For each 
wedge angle a, the arm length L of the barriers was chosen 
such that the length of the longest possible deterministic tra- 
jectory leading through the gap after scattering is equal to 
the persistence path-length L p of the algae (Fig. 3B). We 
then systematically quantified the rectification efficiency by 
scanning the (c?g,ck) parameter space (Fig. 3D). Here, the 
rectification efficiency is defined by the ratio R = (N4) / (Ni), 
where (Ni) denotes the steady-state time- average of the num- 
ber of swimmers in compartment i. Both simulations and 
experiments show maximal rectification for parameter values 
(do = 20 /im, a ~ 30°) that are close to the theoretically pre- 
dicted optimum (Fig. 3D). Thus, although the minimal model 
neglects variability in swimming behavior and hydrodynamic 
effects, it captures the main features of the experiments. 
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Fig. 3. Rectification of algal locomotion in microfluidic ratchets. (A) When con- 
fined in a quasi-2D chamber (height 25 /im), wild-type CC125 perform run-and-turn 
motions, moving ballistically (average speed V ~ 30/im/s) on short time scales 
(< 2 s) and diffusively on larger time scales. (B) Ratchet geometry and schematic 
representation of secondary algal scattering. (C) Rectified steady-state for a chamber 
with four compartments (Movie S3). Scale bar 0.5 mm. (D) In both simulations of 
the minimal model (Materials and Methods) and experiments, the rectification effi- 
ciency R = (N4) / (N±) exhibits a maximum near to the theoretically predicted 
optimal ratchet parameters. Each data point (circles) represents an average over 3 
to 5 different experiments. Rectification efficiencies are linearly interpolated. The 
standard deviation for different experiments is less than 30% of the mean value. 
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Generally, our numerical and experimental results support 
the hypothesis that rectification of microorganisms in environ- 
ments with broken reflection symmetry is a "universal" phe- 
nomenon SSJ Ell [38j [391 HQl ST] |42j |43] . The ratchets mimic 
Maxwell's demon but there is no conflict with thermodynam- 
ics due to the non-equilibrium nature of living systems. 



Conclusions 

Our analysis of individual sperm-surface and alga-surface 
interactions shows that, in contrast to prevailing theoreti- 
cal views, multiple flagellar contact determines the surface 
scattering of these eukaryotic cells. The "puller" -swimmer 
Chlamydomonas reinhardtii exhibits a complex ciliary con- 
tact dynamics over the course of the scattering process, which 
spans a large number of successive swimming strokes and leads 
to an almost complete erasure of memory about the incidence 
angle. Direct contact of the flagella with nearby walls may 
also be the main factor in the surface interactions of other 
motile algae, as the general structure of cilia is highly con- 
served across many eukaryotic species. The flagella- induced 
self-trapping of bull spermatozoa and the Chlamydomonas 
mutant mbol suggests that similar steric mechanisms could 
also be responsible for the long surface-residence times of other 
"pusher" -type microorganisms [HE], possibly even in the case 
of bacteria, where the role of hydrodynamic interactions as 
the main determinant of surface scattering has recently been 
questioned [4J [5j [6] . With regard to future theoretical studies, 
our experimental results anticipate that a detailed qualitative 
and quantitative understanding of microbial surface interac- 
tions will require models that account for the elastic properties 
of eukaryotic cilia and bacterial flagella. 

Finally, as an illustration of how empirically measured 
surface-scattering laws can be exploited to control locomo- 
tion of unicellular green algae in biotechnological applica- 
tions, we have demonstrated robust rectification of random 
swimming for wild-type Chlamydomonas algae. In contrast 
to microfluidic rectification devices for rod-like prokaryotic 
(pusher) swimmers [8], the optimal ratchet geometry for 
algae (pullers) exploits secondary scattering. The proposed 
design can be serialized (see Materials and Methods) and 
parallelized to facilitate large-scale microfluidic implementa- 
tion. Combining the results of the scattering analysis and the 
subsequent rectification study suggests the possibility of inte- 
grating different ratchet geometries to create selection devices 
that sort microorganisms according to structure and dynamics 
of their propulsive appendages. More generally, the present 
investigation implies that suitably designed microstructured 
surfaces can yield new diagnostic tools to quantify function 
and response of eukaryotic cilia (mechano-sensing, desynchro- 
nisation after contact, etc.), which may help to improve our 
understanding of transport processes in the respiratory or re- 
productive systems of higher organisms. 
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Materials and Methods 



Sperm sample preparation. Cryogenically frozen bull spermatozoa were purchased 
from Genus Breeding LTD, UK. For each experiment, a sample of 250/iL was thawed 
in a water bath (37° C) for 15s. The sample was washed three times by centrifuging at 



500$ for 5min and resuspending the pellet in a basic medium containing 72mM KCI, 
160mM sucrose, 2mM Na-pyruvate and 2mM Na-phosphate buffer at pH 7.4 [29l . 

Algal growth. Chlamydomonas reinhardtii strains CC-125 wild-type, CC-2347 
shfl-277, CC-2289 lf3-2, and CC-2679 mbol (The Chlamydomonas Resource Center, 
http : / / WWW . chlamy . org) were grown axenically in Tris-Acetate-Phosphate 
(TAP) medium 22 , on an orbital shaker (200 rpm) in a diurnal growth chamber 
(KBW400, Binder). Daily cycle: 14h at 100/iEm _2 s _1 PAR (Fluora, OSRAM), 
lOh in dark, at 24° C. Cells were harvested during the exponential growth phase. To 
achieve large concentrations of highly motile algae (> 10 6 cells/cm 3 ), cultures were 
centrifuged in Eppendorf tubes at 100g for lOmin and, subsequently, fast-swimmming 
algae were selected according to their ability to swim upwards against gravity. There- 
after, fresh TAP was added to the cell pellet, and the tubes were placed in the diurnal 
chamber for at least one hour to allow cells to recover. 

Microfluidics. Quasi-2D channels were manufactured using standard soft lithogra- 
phy techniques 44 . The master mould was produced from SU8 2015 (MicroChem 
Corp.), spun to a 25/im thick layer and exposed to UV light through a high reso- 
lution mask to obtain the desired structures. The microfluidic chip containing the 
channels was cast from PDMS (Sylgard 184, Dow Corning) and bonded to a PDMS 
covered glass slide after oxygen plasma treatment of the surfaces. Microchannels for 
the sperm experiments were designed such that the sperm cells could be injected in 
a region that was spatially separated from the observation area. In the observation 
area, microchannels (100/im wide) were arranged in a zig-zag pattern with 90° cor- 
ners, where the scattering events were imaged (Fig. 1A,B). The temperature was 
measured by a calibrated thermistor, which was inserted into the PDMS chip 2mm 
away from the observation region. To prevent adhesion of the spermatozoa to the 
walls of the channel, BSA at 5mg/mL was added prior the injection into the cham- 
ber. After injection of the sperm sample into the microchannel, motile spermatozoa 
could escape from the injection site into the observation region. The concentration 
of motile spermatozoa in this region did not exceed 1% volume fraction. In both 
spermatozoa and Chlamydomonas experiments, the microbial solutions were intro- 
duced through inlets that were plugged with unpolymerized PDMS afterwards. This 
procedure prevents fluid flows through the chambers and ensures conservation of the 
total number of cells over the course of the experiment. In the Chlamydomonas 
experiments, the concentration of the algae was kept below 2% volume fraction. 
For the rectification studies, each channel was subdivided into four chambers of size 
2mm X 1mm X0. 025mm, separated by wedge-shaped barriers (Fig. 3). We treated 
PDMS surfaces of the channels prior the experiments with 10% (w/v) Polyethylene 
glycol (m.w. 8000, Sigma) solution in water for 30min to prevent adhesion of the 
algae to the walls, and then flushed them gently with fresh TAP. For the given param- 
eters and uniform initial concentration profile across each chamber, rectified steady 
states were achieved after typically 90min (in the dark). 

Microscopy. To identify the swimming characteristics of individual spermatozoa and 
Chlamydomonas, and their scattering distributions (Fig. 1,2), cell trajectories were 
reconstructed by applying a custom-made particle-tracking-velocimetry algorithm to 
image data taken with a Nikon TE2000-U inverted microscope (lOx objective, lOfps). 
The flagella dynamics close to the boundary (Fig. 1A,B and[2j\,B; Movies Sl,S2 and 
S3, S4, S5 S6) was captured with a Fastcam SA-3 Photron camera (500-2000fps, 
40x/NA 1.3 oil immersion and 60x/NA 1.0 water immersion objectives). For sperm 
white light was used and for Chlamydomonas bright field illumination under red 
light (A >620nm) to minimize phototaxis. Sperm-surface scattering angles were 
determined by tracking the cell body up to a distance of 70/im from the corner. 
For Chlamydomonas, incidence and scattering angles 0i n and ou t were obtained 
by measuring the slope of the trajectory at a distance of 20/im in either direction 
from the scattering point (defined as the location at which the distance from the 
wall becomes minimal). In the rectification experiments, algae concentrations in the 
microfluidic chambers were measured by averaging intensity profiles of image data 
obtained with a confocal scanning microscope Zeiss LSM 700 (5x objective). Here, 
the transmitted light PMT mode was used, while exposing the chambers to laser light 
(639nm) at the lowest intensity to minimize phototactic response of the algae. 

Numerical Simulations. We simulated the dynamics of N = 750 self- 
propelled, non-interacting point particles in a 2D box using MATLAB. Box size 
and ratchet geometries where chosen to match the experimental set-up (Fig. 3). 
To account for the finite radius of Chlamydomonas, a virtual layer of thickness 
a — 5 /im was added to boundaries and obstacles. Particles move ba 1 1 istica I ly 
at a constant speed V until undergoing random turns, or colliding with boundaries 
or obstacles. Initial particle speeds V were sampled from a Gamma-distribution 
r(x; /c, v) with parameters k = 4.2, V — 7.3/im/s, obtained from a best-fit 
to the experimentally measured speed distribution. Run-times between successive 
random turns are sampled from independent exponential distributions with mean T. 
New directions after turn events are sampled uniformly from the unit circle. We 
simulated two types of boundary collision scenarios as simplified approximations 
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to the experimentally observed wild-type scattering behavior: deterministic colli- 
sions with an outgoing angle ou t — 16°, and randomized forward scattering. 
In the latter case, the outgoing angles O ut were sampled from a truncated su- 
perposition of three Gaussians distributions $(x) = y^^ —1 \j<&(x', \±j, (Tj) 
while rejecting angles O ut outside of the interval (0,71"). The distributions pa- 
rameters were chosen as (Ai, A2) = (0.48,0.26) with normalization requir- 
ing A3 = 1 — Ai — A2, and (/ii,/i2,M3) = (0.43,0.26,3.06) and 
(ci, (72, a 3) = (0.58, 0.14, 0.55), obtained from a best-fit to the wild-type 
scattering distribution. In the deterministic case, our simulations confirmed maximal 
rectification for Cz ~ 37°. Simulation results in Fig. 3D show averages from three 
runs with T = 1.5 s for randomized scattering. Generally, numerical results were 
found to be qualitatively robust under moderate variations of the model parameters, 
but the optimal wedge angle is sensitive to changes in the scattering distribution. 

Data Resources and Supplementary Movie Information. A MATLAB script of 
the source code, simulation data and experimental raw data, and additional experi- 
mental movies can be downloaded from: 

http://damtp.cam.ac.uk/ user/gold/datarequests.html 

Scale bars are 20/im in Movies SI and S2; 10/im in Movies S3, S4, S5, and S6; 
0.5mm in Movie S7. 

Serialization (Markov Model). At low-to-intermediate algal volume fractions, as 
realized in our experiments, the dynamics of the Chlamydomonas population on the 
microfluidic chip can be described by a Markov model 7 . This allows to estimate the 



total rectification with an increasing number of compartments. Assuming that the 
chip consists of i = 1, . . . , N identical compartments (Fig. 2D), the time evolu- 
tion of the relative concentration Pi(t) of algae in the ith compartment is governed 
by 



Pi = -kRP! + k L p2, 
Pi = k R pi-! - (k L + k R )pi + k L pi+i, 
Pn — k R p N -i - k L Pn-, 2 < i < N 



[1] 



where k R / R denote the rates for transitions to the left/right neighboring compart- 
ment (pi = dpi/dt and J2iPi — 1)- The stationary distribution {pf}, 
corresponding to the eigenvector of the eigenvalue A = of the transition matrix 
Kij , defined by p% — KijPj> ' s obtained as 



i h i ~ 1 
R 



[2] 



The effective rates k^/n can be estimated by fitting p\ to the experimentally mea- 
sured stationary distribution. Assuming kR > kL, Eq. Q implies that rectifi- 
cation increases with the number of chambers N as p% / Pi = {J^rI^l) N 
Assuming detailed balance and diffusive backward flux suggests that kjifkh — 
PR^dG + ds^/dG, where Pr is probability that an alga is guided through the 
gap after having entered the barrier region. From our experimental data, we estimate 

Pr ~ 0.3. 
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